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1 Introduction 

A key process within astronomy is the exchange of mass, momentum, and 
energy between diffuse plasmas in many types of astronomical sources (in- 
cluding planetary nebulae (PNe), wind-blown bubbles (WBBs), supernova 
remnants (SNRs), starburst superwinds, and the intracluster medium) and 
dense, embedded clouds or clumps (e.g., Fig. ^).. This transfer affects the 
large scale flows of the diffuse plasmas as well as the evolution of the clumps. 
While in much theoretical work this interaction has been ignored, its con- 
sequences can be fundamental, as a growing body of literature now shows. 
Indeed, the standard model of th e interstellar medium is based on such ex- 
changes ( McKee fc Ostrikerll977j) . which occur through, for example, conduc- 



tion, ablation and photoevaporation. The injection and mixing of mass from 
condensations into a surrounding supersoni c medium induces s hocks, increas- 
ing the pressure of the flowing medium fe.g. JPittard et all2 0051. This can lead 
to clump crushing and the reduction of the Jeans mass cau sing star formation 
and i s likely to play a role in sequential star form ation (e.g., Elmegrcen^^J^d^ 
1977), and may allow a starburst to develop l|Hartauist. Dvson k. William; 



1997). Radiative cooling is one way in which a starburst might be regulated, 
as it acts to reduce the pressure of the ambient medium once the mass injec- 
tion rate becomes too high. 

Three lengthscales ch aracterize the entrainmen t of material from a clump 
into a surrounding flow IjHartauist fc Dvsonl 11993). The smallest lengthscale 
is associated with the turbulent boundary layer around the clump. On inter- 
mediate scales the material injected into the flow forms a cometary-like tail, 
such as those seen around clumps in PNe. On the largest scales, the material 
is completely mixed into the flow and becomes indistinguishable from it. Un- 
fortunately, and despite huge effort, the effectiveness of the physical processes 
in controlling the interchange of dense and diffuse material remains uncertain, 
in part because of the complexity of the turbulent boundary layers which ex- 
ist between them. In addition, the microphysics which may drive some global 
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Fig. 1. Numerical simulations which demonstrate the effect of mass- loading on a 
'cluster-wind' from a group of early-type stars. In the left panel the wind expands 
into a smooth medium, while in the right panel the ambient medium is clumpy. 
The time in each simulation is identical. The entrainment of mass from the clumps 
increases the density of the diffuse gas in the bubble interior, reduces its temperature, 
and slows its expansion. 

processes is poorly understood. For example, magnetic reconnection, which 
may be necessary in order for clump and diffuse material to fully mix, is a 
difficult subject. 

The influence of John Dyson in the field of mass-loaded flows cannot be 
overemphasized. He was one of the first to study the process of photoevap- 
oration, he was involved in the development of the widely used analytical 
theory for ablatively-driven mass-loading, and in many subsequent works he 
has investigated the effect of mass-loading on a wide variety of astrophysical 
sources. I am grateful that I have had the opportunity to work with him in 
this field. 

In Sec. |21 1 review our current understanding of mass-injection processes. 
Sec. focuses on intermediate-scale structure, while Sec.0]examines the global 
effect of mass-loading on a flow. Sec. |S] concerns the mass-loading of a vari- 
ety of diffuse sources. For an excellent summary of existing theoretical and 
obs ervational studies on the i nterface between clouds and their surroundings 
see lHartauist fc Dvsonl l|l993j) . 

2 Mass Exchange Processes 

Consider a cloud of radius r c , density p c , and mass M c — 4/37rp c rj?, em- 
bedded in a medium of temperature T, density p, velocity v, and pressure 
P = pkT I [im-ft. Let c c and c be the isothermal sound speed in the cold cloud 
and in the hotter surroundings, respectively, and M. be the Mach number of 
the flow relative to the cloud. Mass can be lost from the cloud and entrained 
into the surroundings through three main mechanisms, as discussed in the 
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following subsections. I describe our current understanding of each process, 
detail analytical estimates of the rate of mass-loss, and highlight current un- 
certainties. The mass-loss rates driven by each process are then compared for 
clumps in a variety of different situations. 

2.1 Hydro dynamic Ablation 

Numerical simulations of the interaction of a supersonic wind or a strong 
shock with a single cloud have been presented many times. The evolution for 
the case of an adiabatic cloud can be broken into 4 consecutive stages: an 
initial transient stage when the shock first strikes the cloud, a compression 
stage, a re-expansion stage, and finally a destruction stage. During the initial 
interaction, a bow shock forms around the cloud, while a shock is driven into 
the cloud with velocity v c w X _1 ^ 2 w s , where x ls the density ratio between the 
cloud and its initial (e.g., pre-shock) surroundings, and v s is the velocity of the 
shock through the ambient medium. The characteristic timescale for the cloud 
to be crushed by the transmitted shock is r cc = r c /v c w X 1 ^ 2 ' r c/ V s- When the 
transmitted shock reaches the back of the cloud, a strong rarefraction is re- 
flected back into the cloud, causing its subsequent re-expansion downstream. 
This is accompanied by a lateral expansion driven by the high pressure in 
the cloud and the lower pressure in the surrounding medium at its sides. The 
cloud is disrupted by the action of both Kelvin-Hcmlholtz (KH) and Raylcigh- 
Taylor (RT) instabilities, with the Richtmyer-Meshkov instability playing a 
minor role unless the surface of the cloud is irregular. Destruction occurs af- 
ter several crushing ti mes, with the cloud material expanding and diffusing 
into the ambient flow ijKlein. McKee fc Colellal ll994fl. In 3-D simulations, in - 
stabilities drive a richer structure llStone fc NormaiJll992tlXu fc Stonelll995ft . 
Recent laser experiments confirm that the the vortex ring which forms at 
the back of t he cloud is broken up by the action of azimuthal bending mode 
instabilities l|Klein et alJl2003|) . In contrast, radiative clouds break up into 
numerous dense cold fragments which survive for many dynamical timescales 
l|Mellema. Kurk & Rottgerindl2002t iFragile et al]l2004h . Self-gravity can be- 
come dynamically important in the dense fragments behind the compression 
shock. External magnetic fields generally increase the compression of the cloud 
and enhance radiativ e cooling, while ma gnetic fields internal to the cloud re- 
sist compression fsee lFragile et al.ll2f)f)5l and references therein). An example 
of a numerical calculation of the time evolution of a cold cloud interacting 
with a supersonic wind is shown in Fig. [21 

An analytical theory for the hydrodynamic ablatio n of material from dens e 
clumps into the surrounding flow was presented by lHartauist et al.l l)l986(l . 
First, consider a clump embedded in a subsonic flow. The magnitude of the 
pressure variations over the surface of the c lump, created by the well-known 
Bernoulli effect, is llLandau fc Lifshit3ll959^ 
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Fig. 2. The destruction of a cold cloud in a supersonic flow by hydrodynamic 
ablation. The time evolution of the logarithm of mass density is shown (in units of 
gcm~ 3 ), with distances given in pc. At the beginning of the simulation the cloud 
center is at z = 20 pc. \ = 500 and M — 3. There is no thermal conduction or 
photoevaporation, magnetic fields, or self-gravity, and the si mulations are performe d 
using 2D axisymmetry. Radiative cooling is included. fFrom lMarcolini et all feoOS), 



\AP\ « R 
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where the stagnation pressure, P s = P + pv 2 . In the small Mach number 
regime, \AP\ as Ai 2 P s . As the flow is fastest at the sides of the clump, the 
pressure is reduced, and the clump expands in directions normal to the flow 
at a speed 

- - ^ hJl+X^M. (2) 



^cxp 



7-1 



For small Mach numbers, v cxp w (7 /2)c c A4 2 ~ c c M. 2 . Mixing between the 
cloud material and the flow occurs within a region of size I which is set by the 
requirement that the rate of mass-loss from the clump, M a b, is comparable 
to the mass- flux of the ambient flow through this region, M s . Since M a b ~ 
M c /t - Af c v oxp /Z, and M s = pvl 2 , I w (M 2 M c C c /pv) 1/3 . In contrast, if the 
flow is supersonic, mixing occurs largely as a result of a low pressure region 
over the reverse face of the clump, 'shadowed' from the flow. Since the mass 
cannot leave the clump faster than its sound speed, v cxp ~ c c , and in this case 
/ « (A/ C c c /H 1/3 - 

The rate of mass-loss from the clump, M a b ~ l 2 piv CKp , where p\ is the char- 
acteristic density of ablated material at distance I. Momentum conservation 
requires that pifexp = pv, so in subsonic flows M a b ~ M 4 ^ 3 (M c c c ) 2 ^ 3 (pv) 1 ^ 3 , 
while in supersonic flow M a b w (M c c c ) 2 ^ 3 (pv) 1 ^ 3 (i.e. independent of M.). 
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These estimates have received some limi ted support from the numerical sim- 
ulations calculated bv lKlein et al . (1994), although the predicted scaling with 
the flow parameters remains to be confirmed. 

An alt ernative approach ba sed on 'mixi ng-length' theory ha s been pre- 
sented bv lCanto fc Ragal l|l99lh (see also lArthur fc Lizanolll997|) . While the 
boundary layer around the cloud is l ikely to be turbulent, eve n if the cloud 
and the surroundings are magnetized ijHartauist fc Dvsonll9 88'l. such theories 
are complicated by the unknown degree to which clump gas and the tenuous 
plasma physically mix, and I do not discuss them further here. 

Finally, it is unclear whether the ablation process by itself can merge 
the stripped material with the global flow in the sense that its temperature, 
velocity, and density approach those of the surrounding tenuous material. It 
may therefore be necessary to invoke another process, such as the transfer of 
heat by thermal conductivity, for the stripped material to acquire the physical 
state of the surrounding medium. Thermal conduction can accomplish this 
phase transition without microscopic mixing, and acceleration to the global 
flow speed is effected by the response of stripped material to pressure gradients 
and viscous coupling, which may arise from a host of mechanisms including 
turbulence. 



2.2 Conductively-Driven Thermal Evaporation 

Cold clouds embedded in a hot medium may also lose mass to their surround- 
ings as hot electrons deposit energy in the surface regions of the clump. This 
process is referred to as thermal evaporation. The rate of mass loss is depen- 
dent on many factors, including the temperature of the hot phase, the clump 
radius, whether the conductivity is saturated, the presence of magnetic fields 
and plasma instabilities, and whether there is a velocity difference between the 
ambi ent medium and the cloud (e.g.. lCowie fc McKedll977t iMcKee fc Cowid 
Nonspherical clumps may be treated in an approxim ate way by adopt- 
ing h alf the largest dimension as the radius of the clump l|Cowie fc Songalial 

US. 

If the mean-free-path for electron-electron collisions, A cc , is approximately 
less than the temperature scale-height, T/|VT|, then the heat flux into 
the cloud, q, is given by the classical theory of thermal conduction i.e. 
q = q c i = -kVT. The mean-free-path A ee = tppjSkT/rrip) 1 / 2 , where the 
electron-electron equipartition time is given by l|Snitzerl ll962') 

_ 3m c 1/2 (fcT) 3 / 2 

47ri/2 nee 4i nA ' [6) 

where In A = 29.7 + In (T/10 6 y/n^) is the Coulomb logarithm and the other 
symbols have their usual meaning. I have implicitly assumed that T c = T. 
The thermal conductivity, n, in a fu lly ionized hydrogen plasma is (see, e.g., 
ISpitzerlll96ilCowie fc McKedll977l) 
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rpb/2 

k = 1.84 x 10~ 5 - — - erg s" 1 K _1 cm" 1 . (4) 
In A 

(the zero current requ irement reduces the ef fective coefficient of conductivity 
by a factor 0.4 - see ISnitzer fc Harmlll953h . The evaporative mass-loss rate 
from a single clump llCowie fc McKee!ll977l) is then 

M con = = 2.75 x 10- u r pc if 2 g s" 1 (5) 

where r pc is the clump radius in parsecs, and T 6 = T/10 6 K. For classical 
evaporation, w = 1. As conductively driven evaporation has a very temper- 
ature sensitive rate, ablation is likely to regulate clump dispersal in lower 
temperature media. 

When A oc > T/|VT|, the classical theory of thermal conduction, which is 
based on a diffusion approximation, may no longer be used. Instead the heat 
flux reaches a limiting value; i.e. it becomes saturated. The approach to sat- 
urated conduction is still partly empirical, and it is c ommon practice to tak e 
a flux-limited form for the heat flux: q sat = 54>pc 3 (jCowie fc McKedll977ft . 
where p and c are in the hot phase, and </) is a parameter of order unity which 
describes the uncertainty in the numerical value of the saturated heat flux. 
Observations of the highly saturated solar wind, laboratory plasm a experi- 
ment s, and Fokker-Planck calc ulations suggest that 0.3 < <fi < 1.1 ((Giulianil 
(1982) conjectured that the effective heat flux, q, can 
be approximated by the harmonic mean of q sa t and q c \, 

111 

-« + — . (6) 

q g sa t q c i 

The resulting heat flux reduces to the smaller of the two conduction forms 
when there is a large disparity between them, and has the convenient property 
of a smooth transition from diffusive to flux-limited transport. The ratio of 
classical to saturated heat flux is 

a = — - K dT (7) 

where the last expression explicitly assumes spherical symmetry i.e. VT = 
dT/dr. The heat flux is then 

k dT 
1 + a dr 

where a is the local stauration parameter. This expression becomes the diffu- 
siv e flux when a <C 1 , and the saturated flux when a ^ 1. 

ICowie fc McKeel il977h also defined a global saturation parameter, 

2kT _( T \ 2 1 
CT ° = 250pc 3 r c ~ V1.54 x 10 7 ) nr pc <j>' ^' 
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where In A = 30 has been assumed. Whereas a measures the saturation lo- 
cally, er measures global scales and allows a quick assessment of the impor- 
tance of saturation effects (<Jq is essentially the local saturation parameter, a, 
evaluated at the ambient conditions with dT/dr — T/r c ). For <7q < 0.03/</>, 
radiative losses quench the evaporati on, and the clump grow s in mass as sur- 
rounding material condenses onto it IjMcKee fc Cowit I II 9771) . This, of course, 
may make the clump gravitationally unstable, and initiate new star formation. 
For 0.03/0 < ao ^ 1, the clump is evaporated at the classical rate. The onset 
of saturation occurs when ao is of order unity, with highly saturated flows 
having cto 3> 1. T he mass-loss ra te in the saturated regime is specified with 
u s» (1 + ero)~ ' 7 (|Giuliani|ll984|) . Since the onset of saturation is dependent 
on the radius of the clump for a specified hot phase, larger clumps will tend 
to evaporate in the classical limit, while the evaporation of mass from smaller 
clumps will tend towards saturation. If the temperature and density of the 
hot phase is evolving (e.g., because it is the interior of a WBB or SNR), the 
radius of clumps which are just at the onset of saturation will als o change. 

Th e dyn amics of the evaporation process have been analyzed bv lMcKee fc Cowie] 
ljl975h and ICowie fc McKeel l)l97/t . and an analogy with ionization fronts 
can be made. For classical evaporation (i.e. 0.03 < ao/4> ~ l)i the veloc- 
ity of the conduction front, -u con d ~ 2a (f>cl/c. If a > 0.25/</>, the con- 
duction front drives a shock into the cloud. Otherwise, the cloud evapo- 
rates subsonically. Conduction fronts in the saturated regime have a velocity 
^cond ~ 1.12(Tq^ 8 c^/c w (for (j> = 1), and always drive a shock into the cold 
cloud ahead of the conduction front itself. 

Time-dependent hydrodynamical simulations of clouds overrun by a strong 
shock and undergoing conductively-d riven thermal evapor ation in the clas- 
sical regime have been calculated bv lOrlando et all l|2005() . Conduction in- 
hibits the growth of Rayleigh- Taylor and Kelvin-Helmholtz instabilities, and 
the fragmentation of the cloud, but heats the evaporated material so that 
it quickly becomes part of the ambient flow. Simulations in the saturated 
regime confirm that the conduction front initially drives a shock into the 
cloud l|Ferrara fc Shchekinovl Il993l iMarcolini et alJ feoOfift . substantially in- 
creasing its density. The resulting evolution then appears sensitive to as- 
sumptions concerning the cooling, with the cloud and the eva poration either 
settling into a quasi-st eady-state llFerrara fc Shchekinovlll9 93'). or displaying 
oscillatory behaviour ijMarcolini et al.ll2005|) . The mass- loss rate during the 
initial phase is comparable to that inferred from Eq. but it is substantially 
reduc ed as the cloud is compressed and decreases in size l)Marcolini et al.l 
2005). Somewhat different mod els includ ing self- gravity have been calculated 
bv lVieser fc Henslerl l)200fjj) and lHensler fc Vieserl l)2002|) . 

Conductively-driven evaporation is perhaps the best-studied of the three 
processes highlighted in this section, but there are uncertainties in many 
physical processes whose influence on M con remains poorly quantified. For 
instance, when ao <; 100, viscous stresses have the potential to increase M con 
significantly, but the exact enhancement depends on the uncertain degree 
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Fig. 3. As Fig. [5] but with heat conduction in the saturated regime (cro = 16). A 
strong shock is driven into the cloud ahead of the conduction front and compresses 
it from all sides, so that its radius has decreased by a factor of 3 by t = 0.3 Myr. 
The subseq uent evolution is also m arkedly different from the purely hydrodynamical 
case. (From iMarcolini et alJl200^ . 



to which they might also saturate IjDraine & Giuliani! 1984). On the other 
hand , if the electron mean free p ath is reduced, the heat flux will be inhib- 
ited. iBandiera fc Chenl l(l994albf l have emphasized that when cro ^ 100, the 
requirement of a zero net current means that hot electrons are stopped by 
electrostatic effects within a thin surface layer. The resulting heat flux is then 
considerably lower than that obtained if the electron s penetrate deep into th e 
cloud and directly heat it through Coulomb collisions l|Balbus fc McKeell982l) . 
The electron mean free path can also be reduced as a re sult of scattering by 
plasm a instabilities, such as the ion-acoustic instability llGaleev fe Natanzori 
1981 . and by whistler waves ijGarv fc Feldmanl 1977t l ^evinson fc Eichlerl 
1992). However, it is difficult to obtain a self-consistent model of these pro- 
cesses, and the presence of a strong magnetic field may suppress such instabil- 
ities. Partial ionization and no n-equilibrium cooling are othe r possibilities for 
suppressing the conductivity IjBohringer fc Harta uist 1987). Magnet ic fields 
may not reduce the conductivity as much as previously thought fe.g.. lBalbusl 
ll986HRosner fc Tuckerll989tlCho et~al"1l2003j) unless the cloud is magnetically 
disconnected from its surroundings. An example of a numerical calculation of 
the time evolution of a conductively-evaporating cold cloud interacting with 
a hot supersonic wind is shown in Fig. |31 



Mass-Loaded Flows 9 



2.3 Photoevaporation 

The fate of a neutral clump exposed to a strong ionizing radiation field has 
been extensively studied over the years. Small, low-mass clumps are instantly 
ionized and rapidly dissipate i.e. are 'zapped'. In contrast, clouds which are 
sufficiently dense and large t rap the ionization front, which becomes D-critical 
jDvsorJll96Sl lBertoldil fl989) and moves into the cloud driving a shock front 
ahead of it. The ionized gas streams away perpendicular to the ionization 
front and expands supersonically into the interclu mp medium, reaching an 
asymptotic Mach number of w 2 (e.g.. lKahnlll969|) . This flow may absorb a 
large part of the incident ionizing flux and appear as a bright rim to the clump. 
The shock driven into the clump is focussed o nto the clump axis, and can 
subst antially increase the initial clump density ijSandford. Whitaker fc Kleml 
GUI). This 'radiation-driven implosion' may make the clump gravitationally 
unstable and lead to new star formation. Otherwise, the pressure overshoot 
causes the clump to reex pand, and it undergoes several radial oscillations 



I Lcflo ch fc, Laz arefi 1994) be fore o btaining an equilibrium, cometary-shaped 
structure l|Bertoldi fc McKeelll990|) . Recombination may occur in the shadow 
of the clump, but this is prevented if there is a diffuse component to the 
radiation field l|Pavlakis et al.ll200l|) . The pressure in the evaporating flow 
declines rapidly, and eventually a termination shock forms. If the surrounding 
medium is supersonic a bow-shock is also formed. Photoevaporated flows also 
occur from th e neutral disks which surround pre-main-se quence stars, known 
as 'proplyds' llO'Dell. Wen fc Hulll993HBallv et al.lll99^ . 

Analytical equations for the mass injecti on rate of the phot oevaporated 
flow a s the clump is destroyed are presented in lBertoldilljl989tl and lMellema et al.l 



(1998). Good agreement with results from numerical simulations is obtained. 
However, it is possible to obtain a simple estimate by setting Af p h = mFA, 
where m is the mass per particle of the neutral material, and F and A are 
respectively the rate per unit time per unit area at which hydrogen ionizing 
photons reach the ionization front and its area. F is reduced by absorptions 
in th e photoevaporating flow, and is approximately given by jMellema et alJ 
1998) 

F « ^ — , (10) 

(1 + a B F r if /3c2) 1/2 

where «b = 2.6 x 10~ 13 cm 3 s -1 is the case B recombination rate for H, m 
is the radius of curvature of the ionization front, q is the isothermal sound 
speed of the ionized gas, and Fq = S /And 2 is the flux delivered by the ionizing 
sources at the position of the clump. Typical assumptions are r;t « r and 
A w 7rr 2 . M p h declines with time as the clump is destroyed and r decreases. 

Photoevaporation may be suppressed if the ram or thermal pressure of 
the surr ounding medium is greater than the pressure of the evaporating flow 
l|Dvsonlll994T) . Density inhomogeneities within the clump may affect the pro- 
cess of photoevaporation, as may the instability of inclined ionization fronts 
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llWilriamsl [ 20021. The role of magnetic fields on the structure of ionization 
fronts jwilharnset al-l EoOOl may also affect M p h- 

2.4 Comparison of Mass-Injection Rates 

While the exact rates of mass- loss by the three processes described in Sec. l2.lt 
12.31 remain somewhat uncertain, an appreciation of their relative importance 
can be obtained by considering several different objects. 

The Helix Nebula (NGC 7293) 

The Helix Nebula is famous for containing many extended cometary tail- 
like structures, which emanate from dense, neutral globules, and point away 
from the central star. The clumps have on av erage the following properties: 
r c w 1CT 3 pc, T c w 10 K, n c « 10 6 cm" 3 ijDvson et alJIlQSit . A typical 



distance of a clump from the central ionizing star is d « 0.1 pc. The star 
has an ionizing photon flux 649 = S/10 49 s» 3.5 x 10~ 4 . The clumps appear 



n 



to be overrun by [OIII] g as with T w 10 4 K 
velocity, v w 17 kms -1 llMeaburn et al-lfeoOSj) . Hence, the Mach number of 
the flow relative to the clumps is about 1.5. Using the analytical equations in 
Secs.EUthrough toE3 we determine that M ab » 1.6 X 10 16 g s _1 and Af ph m 
2.4 x 10 17 g s _1 . The global saturation parameter, ao w 4 x 10~ 7 for </> = 1, 
and gas would like to c ondense onto the clumps at a rate M«2x 10 16 g s _1 
llCowie fc McKedll977ft . though this is likely prevented by the mass-loss that 
occurs through photoevaporation and ablation. The estimated lifetime of such 
clumps is in excess of 4 x 10 4 yr, compatible with the estimated age of the 
nebula. The origin of the knots is discussed in lDvson et al.Nl98<fl . 

The Wolf-Rayet Nebula RCW 58 

RCW58 is a nebula surrounding the Wolf-Rayet star WR40, and which 
has been formed by the current stellar wind sweeping up win d material 
from previous ev olutionary stages, some of which is clumpy ijChul Il982t 
ISmith et alJ [T988) . The clumps are ionized by the central star and typically 
have the following properties: r c ~ 0.1 pc, T c ss 10 4 K, n c w 10 3 cm -3 
( Arthur. Hennev fc Dvsoni riQQfi'l. They are embedded in shocked stellar wind 



material, with n w 1 cm -3 , v « 200 kms -1 . The flow around the clumps 
has M. w 0.6. The H-ionizing flux from the star is ^49 w 2.5 (P. Crowther, 
private communication). For clumps at a distance of 1 pc from the star, 
Co ~ 4, and the conductively-driven mass-loss is mildly saturated with 
M con w 3 x 10 20 gs^ 1 . A higher mass-loss rate is obtained for ablation: 
M a b w 10 21 g s" 1 . Since the cloud is already ionized it does not make sense to 
calculate a photoevaporation rate, but this may once have been the dominant 
process if the cloud was formerly neutral. Therefore, ablation would appear 
to control the rate at which mass is currently stripped from the clump. 
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Within a WBB 

Of course, the mass-loss rate from each process varies within a bubble as the 
density, velocity, temperature, and distance from the ionizing source change. 
This is demonstrated in Fig. where the top panel shows an example of the 
density and temperature structure within a bubble, while the mass-loss rates 
and lifetime of a specific clump are shown in the bottom panel. 





Fig. 4. Top: The internal structure (total number density and temperature) of a 
WBB of age 0.1 Myr expanding into a stationary medium with tih = 10 cm" 3 . The 
central star has M = 10" 6 M Q yr" 1 , Vao = 2000 kms" 1 , and S49 = 1.2. Bottom: 
The mass-loss rate from clumps of radius 0.01 pc, n = 10 r cm" 3 , T — 100 K, 
due to hydrodynamic ablation (ab), photoevaporation (ph), and conductively-driven 
evaporation (con). No clumps are assumed to reside within the central 0.1 pc of 
the bubble. Clumps located within the hypersonic stellar wind are assumed to be 
surrounded by a sheath of hot gas bounded by a bowshock on their windward surface. 
The conductively-driven evaporation is highly saturated for clumps within the hot 
gas of the bubble given the paramaters chosen, but is in the condensational regime 
within the shell of swept-up ambient material. Also shown is an estimate for the 
lifetime of the clumps: t = (M a b + Mean + M p h)/M c ). 
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3 Intermediate-Scale Structure 

While the stripping of mass from clumps has been extensively, though not 
definitively studied, to date there has been very little work on how intermediate- 
scale structures disperse/merge into the background flow. This problem has 
received some a ttention in studies of the interaction between multiple clouds 
and a flow fe.g..| .Tun. Jo nes fc Normanll99fit IPoludnenko. Frank fc B lackman 
l2002USteffen fc L6neJ2004^ . but our understanding of such interactions is still 
developing. A limitation of these works is that the clumps have been modelled 
as single-phase entities, and the density contrast between the clumps and their 
surroundings has, for numerical reasons, typically been ~ 10 2 . The simulated 
clumps then have such short lifetimes that they are unable to significantly 
mass-load the flow. In reality, in many astropysical flows the density contrasts 
are much larger. A different numerical approach which accounts for the much 
longer lifetime s of clumps is to set up sources which continuously inject mass 
into the flow (|Falle et all 12002). Multiple sources act as an efficient barrier 



to the flow if they are sufficiently close together that their combined mass 
injection rate is comparable to or exceeds the mass fl ux of the incident flow 
into the volume that they occupy l|Pittard et al.ll2005l see also Fig. EJ) . In such 
cases, the thermal pressure of the flow is greatly enhanced (at the expense of 
its ram pressure) , and crucially becomes relatively uniform - these conditions 
are exactly those required to increase the probability of cloud collapse and 
new star formation. 

Perhaps one of the best examples of intermediate scale structure is the 
comet -like tail extending from the Galactic Centre source IRS 7, a red s uper- 
giant (|Yusef-Zadeh fe Morrislll99lt fSerabvn. Lacv fe Achtermannlll99l[) . Al- 
though originally interpreted in terms of a global wind- wind collision between 
the slow dense wind from IRS 7 and a Galactic wind, this model had difficulty 
in explaining the length of the tail l|Yusef-Zadeh fc Melialll992]) . However, it 
may be overcome if the wind of IRS 7 is clumpy enough that it becomes semi- 
porous to th e Galactic wind, with mu ch smaller bow-shocks forming around 
each clump ijDvson fe Hartauistlll994|) . 



4 The Global Effect of Mass-Loading on a Flow 

The effect of mass-loading on a global flow can be studied in a one-fluid ap- 
proximation if the following two assumptions are made: i) the clumps are suf- 
ficiently numerous that they can be considered to be continuously distributed; 
ii) the characteristic scale length of injection and mixing is much smaller than 
the dimensions of the global flow (so that the injected mass reaches the gen- 
eral flow velocity and temperature essentially instantaneously) . For the steady 
injection of mass from an ensemble of clumps into the interior of a spherically 
expanding flow, the time-independent continuity and momentum equations 
are then 
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Fig. 5. The interaction of a hypersonic flow with multiple cylindrical mass sources 
when the combined mass injection rate is comparable to the ma ss flux of the inciden t 
flow through the volume that the mass sources occupy ffrom Ipittard et al.ll2005l) . 
The right panel shows an enlargement of the flow around two of the mass sources. 
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where r is the distance from the star, and a is the isothermal sound speed of 
the flow. Q is the mass loaded into the flow per unit time per unit volume. 
The clumps are assumed, on average, to be stationary with respect to the 
global flow, so that there is no net rate of momentum injection. 

The continuity and momentum equation may be combined to give 



dr 



Q f 2 , 2\ , 2a2v 

— {v + a ) H . 

p r 



(13) 



If r is large (i.e. the flow is plane-parallel) and if v > a (i.e. M > 1), 
then v 2 > a 2 , dv/dr < 0, and the flow decelerates. On the other hand, 
if r is large and v < a (i.e. M < 1), then v 2 < a 2 , dv/dr > 0, and 
the flow accelerates. A key feature of mass- l oading is that it tend s to drive 
the ambient flow to Mach number unity l Hartauist et alJ flQSfi'l. For an 
expanding, tim e-dependent flow, mass-loading m ay drive the Mach num- 
ber to 0.6-0.7 ijArthur. Dvson k Hartauistl Il993|) . Another feature is that 
if mass-loading is large, a global reverse shock will be greatly weakened 
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llArthur. Dvson & Hartauistll994 


Williams. Hartauist & DvsorJl 99511 Williams. Dvson & Hartauistl 


Il999l IPittard. Hartauist & DvsoiJ 


2001b). In addition, a two-fluid study has 



reveal ed that shocks which ov errun clumpy media are weakened and broad- 
ened {Williams fc Dvsonll2002]) . This has significant implications for the sur- 
vival of the clumps. Finally, the injection of mass into a radiative medium has 
a stabilizing influence against isobaric and isentropi c perturbations, and can 
suppress the development of the thermal instability IjPittard et al. 2003) . 



5 Sources with Large-Scale Mass-Loaded Flows 

Theoretical studies of the effect of mass-loading on a flow fall into two cate- 
gories. If the flow is treated as a single fluid and has a very simple geometry, 
it is often possible to obtain a similarity solution. Alternatively, one can use 
a hydrodynamical code to model and find a wider variety of solutions. This 
allows much greater flexibility, and the ability to incorporate a larger num- 
ber of physical processes. In principle, a two-fluid approach can be used, in 
which a time-dependent spectrum of clumps is separately tracked. Of course, 
an unavoidable drawback of such computations is the loss of generality, as, 
for instance, the clump spectrum inside a starburst superwind is likely to be 
somewhat different to that in our local ISM. 

In the following subsections I discuss mass-loading in a variety of as- 
trophysical settings. The properties of winds mass-loaded by material from 
chimps (or stellar sources) has al so received theoretical attention in the con- 
text of ultracompact H II regio ns lIDvson. Williams & Rcd manll995lh Herbig- 
Haro objects (see discussion inFH artquis t fc Dvsonlll988lll993|) . and globular 
cluster winds {Durisen fc Burndll98lj) . The mass-loa ding of accretion flows in 
stell ar clusters and AGNs has also been s tudied bv lDavid fc Durisenl |l989) 
and lToniazzo. Hartauist fc Durisenl l 200ll) . 



5.1 WBBs, PNe, and Superbubbles 

The dramatic and beautiful impact of stellar winds on their environment has 
long been recognized (e.g., see the chapter by Jane Arthur). The first simi- 
larity solutions of this inte raction, which assumed a smo oth wind and envi- 
ronment, were obtained bv lDvson fc de Vriesl l|l972l) and lDvsonl jl973h . The 
shocked ambient gas was fo und to cool very rapidly (as subsequently shown 
in a numerical simulation bv lFalleill975l . and the hot shocked wind was found 
to fill most of the bubble vo lume. The effect of mass-loading was first consid- 
ered bv lWeaver et all l|l977h . who obtained an approximate similarity solution 
for a WBB expanding into a constant density medium by assuming an iso- 
baric shocked wind region. In this work, mass-loading of the interior of the 
WBB is assumed to occur through the conductively-driven evaporation of the 
cool swept-up shell. The mass transferred into the bubble interior m a y eas- 
ily dominate the total mass within the bubble. lHanami fc Sakashital |l987h 
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M1-67 




Fig. 6. The Ml-67 nebula around WR124. This is the youngest known WR neb- 
ula, and unlike older nebulae there is no obvious global shell surrounding a central 
cavity. Instead, the interaction between the current wind and its LBV progenitor 
starts close to the cent ral star and extends to the outer boundary (reproduced from 
iGrosdidier et al.lll99gft . 

used this same approach to generalize the I Weaver et al.l l(l977h model to an 
ambient density with a power-law slope, as well as considering the conduc- 
tive evaporat i on of embedded clumps. An analytical solution based on the 
Weaver et al.l lll977l) model for W BBs in a density gradient was obtained by 
Garcia-Segura fc Mac Low! l)l995(l . Similarity solutions of WBBs mass- loaded 
by embedded clumps and obtained without the assumption of a constant pres- 
sure shocked wind we re described bv lPittard. Dvson fc Hartauistl il2001al) and 
IPittard et all I^OOlbj) . A similarity solution analogous to the I Weaver et alJ 
lll977l) model, but t aking into account the tim e-depe ndence of the stellar 
wind, was obtained bv lZhekov fc Perinottol ||199(t| ) . whilelZhe koy fc Mvasnikovl 
(1998) derived a similarity solution (based on the lWeaver et al.l l|l977h model) 
with reduced thermal conduction. 

Some of the best studied WBBs are around WR stars. Many of these 
display a clumpy morphology, and mass-loading is likely to be an important 
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proces s. Evidence for mass-loading in WBBs was obtained by ISmith et al.l 
l|l984|) . who observed a correlation between the velocity and ionization poten- 
tial of ultraviolet absorption features towards the central star of the RCW 58 
WBB. The slope of this relationship effectively measures the ratio of the mass 
flux to the pressure in the line- forming region, i.e. pv/P, and the value of 
this ratio implies that the rate of mass loading into the bubble is 40-50 ti mes 
greater than the mass- loss rate of the WR star ((Hartauist et al.lll98f)() . A 
standard bubble without mas s-loading is unable to explain the wide range of 
observed velocities, while the IWeaver et al ] lll977l) model predicts a correla- 
tion of the opposite sense, irrespectiv e of the assumed ambient density gra- 
dient (see lHanami fe Sakashitalll987f) . In contrast, similarity solutions with 
mass loading from embedd e d clumps are able to r eproduce this correlation 
l(Hanami fc Sakashitalli987t IPittard et aHl2001albh as ca n hydrodynamical 
simulations of the same process llArthur et al.lll993l Il996|) . Further evidence 
for mass-loading in wind-blown bubbles is provided by spectropscopic data of 
core-halo PNe, which indicate halo tem peratures in excess of t hat wh i ch can 
be obtained by photoio nziation alone IjMeaburn et alJ Il99ll) . IDvsotiI ill 9921) 
and I Arthur et all l|l994h have shown that this is consistent with a transonic 
wind leaving a mass-loaded core region and shocking against clumps in the 
halo region. 

Without some form of mass-loading, the hot gas in WBBs, PNe, and super- 
bubbles may be too rarefied to produce observable X-ray emission. Irrespec- 
tive of the entrainment process, mass-loading increases the interior density of 
the bubble while simultaneous reducing its temperature (the same thermal 
energy must now be shared between more particles), and as a consequence, 
the X-ray emission both softens and dramatically increases. We are therefore 
fortunate that the last three decades have been a golden age of X-ray astron- 
omy, with the operation of successive facilities with vastly improved sensi- 
tivity, spatial- and energy-resolution. The observed temperatures in WBBs, 
PNe, and quiescent (i.e. those without recent SNe) superbubbles are all rea- 
sonably soft and therefore indicative of mass-loadin g (though oth e r exp lana- 
tions are also possible). However applications of the lWeaver et alJ l)l977|) and 
ICarcia-Segura fc Mac Low! l|l995l) models predict far too much X-ray emission 
(sometimes by a factor of 100), so at the very least the evaporation of mass 
from a dense swept-up shell must occur at a much reduced rate. 

Only two WBBs hav e been detected t o date at X-ray energies, NGC 6888 
l)Wriggelll999() and S308 l|Chu et al Jl2003^ . both of which are limb -brightened. 
At firs t, thi s morphology seems to agree with th e predictions of the lWeaver et all 
l|l977t) and lOarcia-Segura fc Mac Lowl 1^)95) models since the emissivity is 
highest where the density is highest and the temperature lowest. However, 
the soft emission from the limb is easily absorbed by t he ISM, with the re - 
sult that such models predict a centre-filled appearance ijWrigge et al.ll2005|) . 
Again, a reduction in the conductivity can improve the level of agreement. It 
remains to be seen whether bubbles with distributed mass-loading from em- 
bedded clumps can better reproduce the observed limb-brightening, though 
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this process has some support from hydrodynamical simulations which often 
show instabilities breaking dense clumps off cold swept-up shells and sub- 
sequently becoming entrained in the hot bubble interior (e.g., see Fig. 6 in 
the chapter by Jane Arthur). An alternative explanation is that the X-ray 
emission from NGC 6888 and S308 arises from the shocked RSG wind, and is 
curren tly enhanced in these objects b y a collision between the WR and RSG 
shells ()Frever. Hensler fc Yorkell2006l) . However, this requires that the RSG 
wind is fast (~ 1 00 kms" 1 ), and speeds of this order are far from certain 
(see discussion in lGarcfa-Segura. Langer fc Mac Lowlll996|) . Finally, there is 
a vast amount of evidence that the stellar winds themselves are clumpy, which 
has further implications for their interaction (see, e.g., Fig.^J). 

Discerning the morphology of PNe is more difficult, as their angular 
size is smaller. Only Mz 3 and NGC 6543 are adequatel y resolved , and, 
like the WBBs aro und massive stars, are limb-brightened ijChu et alJl200li 
iKastner et aLlEool^ . Unlike the WBBs around massive stars, internal absorp- 
tion within the nebula is likely to be important, while the action of collimated 
outfl ows and heat condu ction may also modify the observed X-ray morphol- 
ogy ijKastner et alJl2002f) . The duration for the presence of hot gas appears 
to be short, as only young PNe show diffuse X-ray emission. 

On larger scales, the combined action of stellar winds and SNe from mas- 
sive stars in OB associations sweep up the ambient ISM to form superbubbles. 
The physical structure of a superbubble is expected to be similar to that of 
a WBB formed around an isolated massive star. Superbubbles containing 
recent SNe have X- ray luminosities which exceed the predictions from the 
Wea ver et al 1 lll977l) model, but when in a quiescent state (i.e. without recent 
SN blasts) the X-ray luminosities ar e an or der of magnitude lower than ex- 
pected from the IWeaver et al ] lll977l) model llChu et aLlllfffll . In contrast to 
WBBs and PNe, quiescent superbubbles have a center-filled X-ray morphol- 
ogy, with brighter X-ray emission near the central star cluster, and hotter 
gas (T ~ 10 7 K). The diffuse X-ray emission from hig h-mass star-forming re - 
gions has also been studied on smaller scales (see, e.g.. lTownslev et al"1 l2003). 
In M 17, it appears that the hot plasma is mass-loaded by a factor of ~ 10 
(Dunne < EHIC2003). However, other processes, such as particle acceleration, 
and different electron and ion temperatures may also be important. 

An interesting finding for clump-embedded, conductively-driven mass- 
loading is the occurence of a negative feed-back mechanism, whi ch sets a max- 
imum limit to the amount that a bubble can be mass-loaded ( Pittard et al.l 
l2001afl . This arises as a consequence of the dependence of the evaporation 
rate on temperature and the lowering of temperature by mass loading. Since 
mass-loading to primary stellar wind mass ratios approaching those inferred 
for RCW 58 could not be obtained, a key conclusion is that ablation from em- 
bedded clumps must be the dominant driver of mass-injection in RCW 58 in 
particular, and perhaps also more generally (this finding is in agreement with 
the estimate that M a b > M con in Sec. 12.4(1 . However, the lower mass-loading 
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to primary stellar wind mass ratios deduced for S 308 and M 17 probably do 
not rule out conductively-driven mass-loading in these sources. 

5.2 SNRs and Starburst Superwinds 

The fact that the ISM is known to be multi-phase means that SNRs are un- 
doubtedly mass-loaded, and there is clear observat i onal evidence of engulfed 
clump s within the SNR N63A (|Chu et alJ Il999t IWarren. Hughes fc Slanel 
2003). Young SNRs first interact with the circumstellar material ejected by 
their progenitor. The youngest SNR in the Galaxy, Cas A, contains bright, 
slow-moving knots of gas called quasi-stationary flocculi (QSFs), which have 
been suggested to arise from the circumstellar bubble of the WN progenitor 
llChevalier fc Kirshnerlll97ah . 

Similarity solu tions for SNRs ma s s-load e d from embedded clouds have 
been obtained bv iMcKee fc Ostrikerl l)l977|) . IChieze fc Lazarefj lll98ll). and 
IWhite fc Lond l|l99lj) for conductivelv-driven evaporation, and bv lDvson fc Hartauistl 
l)l987[) for ablation-driven injection. A small number of papers based on nu- 
meric al simulations of mass-loaded supe rnova remnants also exist in the liter- 
ature. ICowie. McKee fc Ostrikerl l(l98l|) included the dynamics of the clumps 
and found that warm clumps are swept towards the shock front and are rapidly 
destro yed, while cold clumps ar e more evenly distributed and have longer life- 
times. j^^Su^^^enne^ {L996) studied the effects of mass loading by hydro- 
dynamic ablation on supernova remnants evolving inside cavities evacuated 
by the stellar winds of the progenitor stars. 

When SNRs overlap in regions with vigorous star formation, they may cre- 
ate highly pressurized superbubbles that burst out into intergalactic space. In 
the standard picture of such starburst superwinds, t he wind remains subsoni c 
until it reaches the boundary of the starburst region llChevalier fc ClegdlQSffl . 
However, the predicted X-ray luminosity of the thermalized SN and stellar 
wind ejecta is lower t han observed unless the superwind is heavily mass-loaded 
l|Suchkov et al.lll99(T|) . with the rate of mass-injection from t he destruction of 
cloud s several times larger than that due to SNe (see also iHartauist et alJ 
19971) . Recently, high-spatial-resolution observations have provided a much 
clearer view of the mass-loading process in superwinds. Cool Ha emitting fil- 
aments and clumps are seen embedded within the superwind, and the soft 
X-ray emissio n appears to be associ ated with hot gas interacting with these 
structures (see lStrickland et all20o4 and references therein) . Theoretical sup- 
port comes from hydrodynamic models which show that the superwind sweeps 
up and inco rporates large masses of ma terial from within the galactic disk as 
it develops (|Strickland fc Stevensll2000|) . 

Since superwinds are driven by overlapping SNRs, and the range and 
radiative energy losses of a remnant are affected by mass loading, it is 
desirable to have approximations which describe remnant evolution and 
range in clumpy media. The first steps towards this goal were taken by 
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iDvson. Arthur fc Hartauistl l|2002fl and IPittard et all l|2003l) . where hydrody- 
namical simulations of SNRs undergoing mass-loading driven either by abla- 
tion or conduction were calculated. Significant differences between the evolu- 
tion of the SNRs were discovered, due to the way in which conductive mass 
loading is extinguished at fairly early times, once the interior temperature 
of the remnant falls below ~ 10 7 K. At late times, remnants that ablatively 
mass load are dominated by loaded mass and thermal energy, while those that 
conductively mass load are dominated by swept-up mass and kinetic energy. 
Th ese works may ultimately be used in superwind models which are akin to 
the iMcKee fc Ostrikeil ljl977j) model of the interstellar medium. 



5.3 AGN-SNR Interaction 



Many theoretical explanations have been propo sed for the origin of t he broad 
emission line regions (BELR) in AGNs (see, e.g.. |Pittard et alJ2003bl and ref- 
erences therein), i ncluding the interaction of an AGN wind with supernovae 
and star clusters l|Perrv fc Dvsonlll985l IWilliams fc Perrvlll99^ . It is now 
clear that the bulk of the broad-line emission arises from an accretion-disk 
wind ijProga. Stone fc Kallmanl 1200(1 see also the chapter by Stuart Lums- 
den), but interactions between an AGN wind and SNRs may make a non- 
negligible contribution to the emission of high i onization line s. These inter- 
actions will also mass-load the AGN wind (e.g. Smith 1996), and possibly 
could be used as a diagnostic of it. Hydrodynami cal simulations of the early 
stages of this interac tion have been presented bv IPittard et alJ l)2001cf) and 
IPittard et al. (20023)- The strong radiation field means that cool post-shock 
gas forms only for a very limited time, before being heated back up to the 
Compton temperature. 



5.4 Intracluster Gas 



The intracluster medium consists of hot, subsonic gas which is bound within 
the gravitational potential of the cluster. It was suggested almost 3 decades 
ago that the density of gas within the central regions is high enough to permit 
significant cooling within a Hubble time This energy loss causes outer gas to 
flow subsonically in to the centre in order to maintain hydrostatic equilibrium, 
a process referred to as a 'cooling flow'. Recent X-ray observations appear to 
show a systematic deficit of low temperature emission in comparison to the 
standard isobaric cooling-flow model, which has led to the questioning of this 
model. Mass de position rates significa ntly smaller than expected are also in- 
ferred (see, e.g.. IPeterson et aTll2003l and references therein). Attention has 
now focussed on the possibility that the gas is prevented from cooling by some 
compensating form of energy injection. For instance, the frequent occurence of 
X-ray cavities coincident with radio lobes around the central dominant galaxy 
demonstrates that AGNs can significantly influence the X-ray morphology of 
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Fig. 7. An Ha image showing the filamentary structure around NGC 1275 (credit 
C. Conselice (Caltech), WIYN, AURA, NOAO, NSF). 

the hot gas (e.g- lFabian et alJl2002|) . However, the timescale for energy trans- 
fer from the relativistic plasma to the surrounding gas is poorly known, while 
the effectiveness of heat conduction at preventing cooling is enthusiastically 
debated. 

An alternative interpretation is that the cooling plasma radiates its en- 
ergy in the UV/optical bands rather than at X-ray wavelengths. This can 
occur if there is mi xing between the cooling plasma and colder material (e.g., 
iFabian et all 12002ft . and may be reprocessed into the IR if the gas is dusty. 
Luminous optical /UV nebulosity is common in cluster cooling flows (e.g., 
iHeckman et al.l Il989^ . and is particularly widespread a round NGC 1275 in 
the Perseus cluster l|Conselice. Gallagher fc Wvsell200l|) (see Fig. 0). While 
present information is rather sparse, and some of this emission will be pow- 
ered by star formation, the total submillimetrc to UV e mission appears to b e 
sufficient to account for the missing soft X-ray emission l)Fabian et alJ F2002). 
However, direct evidence for the UV e mission resulting fr om this mixing is 
rather tentative in curr ent FUSE data l|Oegerle et al.ll200l|) . 

If cold clouds are embedded within the hot gas, they could be important for 
global mass injection into the hot surroundings. Indeed, since the intracluster 
gas is observed to be enriched with metals, material must either be injected by 
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the host galaxies, through fountain flows or superwinds, or ablatively stripped 
from them. Mass-loading usually reduces the mean temperature of a flow, as 
the thermal energy is shared between more particles, but in clusters this sit- 
uation is modified by the gravitational forces on the flow and the significant 
gravitational potential energy that the entrained material possesses. A pre- 
liminary study has shown that under such conditions it is pos sible to increase 
the radiated power at high temperatures IjPittard et alJl20f)4|) . Hence the dif- 
ferential luminosity distribution can have a positive slope with T, as required 
for agreement with the latest observations of clusters. In addition, ripples 
found in the X-ray emission map of the Perseus cluster l|Fabian et aLl boOS) 
resemble the structure of some mass-loaded accretion flows. The most likely 
explanation for their origin may be time variations in the source of outflowing 
material, but this interaction clearly needs further investigation. 

6 Summary 

The mass-loading of flows is now recognized as a fundamental process in 
astronomy, yet much work remains. We still do not have a good understanding 
of ablatively-driven evaporation or conductively-driven thermal evaporation 
and the ability of plasma instabilities to quench it. The rapidity with which 
cold material is mixed into a hotter flow is also unknown. The global properties 
of a variety of mass-loaded flows have now been studied, but in only a few 
instances have specific objects been modelled and a direct comparison with 
observations made. 
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